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We report the first observation of Xc2 production in B-meson decays. We find an inclusive 
B -> xciX branching fraction of (1.80±q ± °- 26 ) x 10_3 - The data set > collected with the Belle 
detector at the KEKB e + e~ collider, consists of 31.9 million BB events. We also present branching 
fractions and momentum spectra for both Xci and Xc2 production. 

PACS numbers: 13.25.Hw,14.40.Gx,14.40.Nd 



Although the theory for weak decays of b quarks is 
formulated in terms of quark processes, experiments are 
done with B hadrons. The application of quantities cal- 
culated at the quark level to the physically realizable 
hadrons usually requires theoretical assumptions and ap- 
proximations. One widely used approximation is "fac- 
torization," where it is assumed that the participating 
quarks form hadrons with no subsequent transfer of quan- 
tum numbers between them Q . Since this assumption is 
widely used, it is important that the range of its validity 
is carefully tested. 

In the factorization limit, decays of the type B — > XcqX 
and XciX are not allowed by angular momentum and 
vector-current conservation [0. These decays can occur 



if there is a (factorization- violating) exchange of soft glu- 
ons between the quark pairs prior to hadron formation. 
Belle has recently reported the observation of the de- 
cay B~ — ► XcqK~ with a decay branching fraction that 
is comparable to that for the factorization-allowed de- 
cay B~ -> J/ipK~ f|. The CLEO collaboration has 
published a 95% CL upper limit on the inclusive decay 
B -> Xc2X of 2.0 x 1(T 3 f§. 

In this paper we report evidence for the inclusive decay 
B — > Xc2X from an analysis of 31.9 million BB events 
produced in a 29.4 fb _1 data sample taken at the T(4S) 
resonance with the Belle detector at the KEKB asymmet- 
ric e + e~ collider. An additional 3.0 fb" 1 sample taken 
at a center-of-mass energy 60 MeV below the T(45) is 
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used to study backgrounds from non-resonant (contin- 
uum) processes. 

The Belle detector consists of a three-layer silicon ver- 
tex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov coun- 
ters (ACC), time-of-flight scintillation counters (TOF), 
a CsI(Tl) crystal electromagnetic calorimeter (ECL), a 
1.5 T superconducting solenoid coil and an instrumented 
iron- flux return for muon and Kl detection (KLM). The 
detector is described in detail elsewhere || . 

Events with candidate B mesons are selected by first 
applying general hadronic event criteria. These include 
the requirement of at least three charged tracks, an 
event vertex consistent with the interaction point, recon- 
structed center-of-mass (CM) energy greater than 0.2^/s, 
a longitudinal component of reconstructed CM momen- 
tum less than 0.5y/s/c, and a total ECL energy between 
0.1-y/i and 0.8y/s with at least two energy clusters. To 
suppress continuum backgrounds we also require the ra- 
tio of the second to zeroth Fox- Wolfram moments to be 
less than 0.5 jfj). 

We reconstruct Xd an d Xc2 via the decays to J/ipj, 
J/ip — ■» l + l~ . Both leptons are required to be loosely 
identified as leptons. Electrons are identified using 
a combination of drift chamber dE/dx measurements, 
aerogel response, and electromagnetic shower position, 
shape and energy. Muons are identified with KLM hit 
positions and penetration depth. In order to recover di- 
electron events where one or both electrons have radi- 
ated a photon (final state radiation or bremsstrahlung) , 
we include the four-momentum of every photon detected 
within 0.05 radians of the original e + or direction in 
the invariant mass calculation. The J/ip — ► /i + /i~(e + e~) 
candidate invariant mass is required to be between 
-25(-40) MeV/c 2 and +25 MeV/c 2 of the known J/ip 
mass, with an expected resolution of 9.6(10.8) MeV/c 2 
for dimuon(dielectron) J/ip's. The larger range for di- 
electron candidates is to include candidates that fall in 
the radiative tail, even after the photon correction. 

To reduce combinatoric background, we veto gamma 
candidates that form a good ir° candidate with any other 
photon candidate of energy greater than 60 MeV in the 
event. A good candidate n° is defined by a % 2 of less 
than 10 after a mass-constrained kinematic fit. We then 
make a histogram of the mass difference between the Xc 
and the J/ip candidates; this nearly eliminates the effect 
of the J/ip measurement error. The error on the mass 
difference is dominated by the photon energy resolution. 
The momentum of the Xc candidate in the CM refer- 
ence frame is required to be less than 1.7 GeV/c (the 
kinematic limit for a Xc coming from a B meson); this 
requirement was not used in the determination of the Xc 
momentum spectra. 

In Fig. [I] a clear Xc2 peak can be seen next to a larger 
Xci peak. In order to determine the yield we fit the dis- 
tribution to two Crystal Ball line shapes H and a third- 
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FIG. 1: The \c — J/ip mass difference distribution for candi- 
date events. The widths reported correspond to the "width" 
parameter of the Crystal Ball function [J7J . 

order Chebyshev polynomial for the background. The 
Crystal Ball function allows for a "tail" in the line shape 
that is due to photon shower leakage in the ECL. 

In this fit (the "standard" fit), the signal line shapes 
(i.e. the widths, means, and tail parameters) are allowed 
to float with the following constraints: the difference be- 
tween the means is fixed to the known Xci — Xc2 mass dif- 
ference; the Xc2 width is fixed to 1.1 times the Xd width, 
to take into account the Monte Carlo expected ratio of 
the widths, which is consistent with a higher average Xc2 
photon energy; and the tail parameters are fixed to be 
the same. The background shape is fixed by fitting to 
the regions outside the signal region from 0.35 to 0.50 
GeV/c 2 . 

The signal shape was compared with predictions from 
an inclusive B — > XciX and Xc2X full Monte Carlo sim- 
ulation. The signal widths in data are larger. In a study 
of D*° — > £>°7, we find that the calorimeter response for 
a single photon is broader in the data than in the Monte 
Carlo: for the Xc — J/*P mass difference, we expect the 
width to be increased by a factor of 1.3. For Xci the 
Monte Carlo width is 7.0 ± 0.2 MeV/c 2 , the corrected 
width is 9.1 ± 0.3 MeV/c 2 , and the measured width is 
10.0 ± 0.6 MeV/c 2 . We consider the variation in signal 
yields for various fitting scenarios in determining the sys- 
tematic error due to fitting. 

The background shape was checked against a full 
Monte Carlo simulation that included the appropriate 
amounts of BB and non-resonant events. The Monte 
Carlo and data background shapes are in good agree- 
ment and their normalizations agree within 3%. 

We find a yield of 2529 ± 127 events in the Xci peak 
and 611 ± 76 events in the Xc2 peak, where the error is 
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statistical only. 

Several sources of background production were 
checked. Two-photon processes produce Xc2 @< To es- 
timate the contribution to the Xc2 signal from events of 
this type we looked at the equivalent of 560 fb _1 of Monte 
Carlo data. From this sample we estimate a background 
contribution of 1.9 events. We also checked the 3.0 fb _1 
continuum data sample for \c production. We expect a 
small number of events from feed down from continuum 
ip(2Sys and possible direct Xc production. From the fit 
we find 14.0 ± 6.4 events in the Xci region and 0.4 ± 5.7 
events in the Xc2 region. Expected contributions of feed 
down from continuum iJj(2S) production are 0.5 events 
for Xci and 0.2 events for \ci [0, and, hence, consis- 
tent with the above measurements. For the Xc2 case, we 
follow the prescription of Feldman and Cousins and find 
the 68.27% confidence interval for the event yield to be 
[0.0,6.1] (ll|. We scale the continuum yields by the ra- 
tio of on- and off-resonant luminosities, corrected for the 
difference in continuum cross section due to the slight dif- 
ference in beam energies. The scaled Xci and Xc2 contin- 
uum yields are subtracted from the on-resonance yields. 
We use the Feldman- Cousins confidence limits in deter- 
mining the statistical error for Xc2 after the subtraction. 

To convert yields to branching fractions we deter- 
mine the reconstruction efficiency with a full inclusive 
B — > XciX and Xc2X Monte Carlo. We find the efficien- 
cies for reconstruction to be 32.0±0.5% and 33.1±0.9%, 
respectively. The Xc momentum spectra of the Monte 
Carlo are similar to those measured in data. The effi- 
ciencies are uniform over the allowed Xci , Xc2 momentum 
range. 

We use the 2001 Particle Data Group Q values for 
daughter branching fractions B(J/ip — > l + l~) = 0.118 ± 
0.002, B(Xd J/^l) = 0.273 ± 0.016, and B( X c2 -> 
J/ifrf) = 0.135 ± 0.011. The inclusive B -> Xc X branch- 
ing fractions are found to be: B(B — > Xc%X) = (3.63 ± 
0.22) x 10~ 3 , and B(B -► x c2 A) = (1.80±gJ|) x 10~ 3 . 
These numbers are summarized in Table |[ 

Some of the B — > Xc decays result from "feed down" 
from the ^(2S I ); these are not forbidden by factorization. 
In order to determine the rate for direct decays to the Xc 
states, the i/j(2S) contribution must be subtracted. This 
feed down is estimated using the Particle Data Group 
B — > ip(2S)X and ip(2S) — > Xci branching fractions. 
After correcting for feed down we find: B(B — ► XciX) — 
(3.32 ± 0.22) x 10~ 3 , and B(B -> X c2X) = (1.53±g;§|) x 
10~ 3 . 

Significant sources of systematic error are in the effi- 
ciencies for lepton identification (2% per lepton track), 
tracking (2% per track), photon detection (2%), as well 
as daughter branching fractions (6% for Xci, 8% for Xc2), 
and fitting systematics (4% for xa, 10% for Xc2)- The 
systematic errors are summarized in Table H. 

The fit for the Xci and Xc2 yields is sensitive to the 
signal and background shapes. We estimate the error 



associated with the fit by performing the fit in a vari- 
ety of ways including: fixing the signal means, widths, 
and tail shapes to Monte Carlo values (with the widths 
multiplied by a scaling factor and separately by adding a 
random number from a Gaussian distribution generated 
to yield the desired width increase); allowing the means 
to float, with the widths and tail shape fixed; allowing 
the means and widths to float, with the tail shape fixed; 
and allowing all parameters to float. In all cases, when a 
parameter is allowed to float, the Xci and Xc2 hue shapes 
are constrained appropriately as with the standard fit. 
Two methods of fitting the backgrounds are also used: 
fixing the background with the sidebands (as with the 
standard fit) and allowing the background shape to float 
freely. The one combination that is not used is to fit with 
a free tail shape and a free background shape as there can 
be a trade off between the background area and tail area 
in the fit. 

In addition to the above fits, we confirmed that a third- 
order polynomial is sufficient to fit the background by 
performing a fit to the background Monte Carlo; adding 
additional terms did not improve the confidence levels of 
the fits. The fitting systematic error is assigned from the 
largest variation between the fits described above and our 
standard fit. 

The Xc momentum spectra are interesting as they can 
give clues to the production mechanisms. The high 
momentum end is dominated by two-body decays to 
Xci(Xc2)K and Xci(Xc2)K* while the low end may be 
from higher mass K* resonances, multi-body decays or 
feed down from ip(2S). To determine the momentum 
spectra, we divide the data into sets based on the mo- 
mentum of the Xc candidate. We then fit each distribu- 
tion for the Xd and Xc2 yields, which are converted into 
differential branching fractions, corrected bin- by-bin for 
the detector efficiency. The resulting momentum spectra, 
shown in Fig. ^, are broad indicating that a large compo- 
nent of either multi-body decays or higher K* resonances 
is present. The shaded histogram in Fig. || shows the 
Xc2 momentum distribution for Monte Carlo-simulated 
B — > Xc2K decays, which indicates that almost all Xc2's 
from these decays have momenta between 1.2 and 1.6 
GeV/c. After doing a fit of this Monte Carlo histogram 
to the data histogram we find an upper limit at the 90% 
confidence level of 5.0 x 10 -4 for the B — > Xc2K branch- 
ing fraction. The shaded area in Fig. || corresponds to 
this upper limit. A more detailed analysis of this decay 
is forthcoming. 

In summary, we report the first statistically signifi- 
cant observation of Xc2 production in £?-meson decays. 
The B — ► Xc\X and B — ► Xc2X branching fractions 
are measured to be (3.63 ± 0.22 ± 0.34) x 10~ 3 and 
(1.80i^28 ± °- 26 ) x icr3 > respectively, where the first 
error is statistical and the second systematic. After sub- 
traction for feed down from ip(2S), we find the direct 
branching fractions to be (3.32 ± 0.22 ± 0.34) x 10~ 3 and 
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TABLE I: Yields and branching fractions. Errors are statistical only. 
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TABLE II: Systematic Errors. 
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FIG. 2: Branching fractions for B — * XaX and B — > x^X 
as a function of \c momentum in the e + e~ center-of-mass 
frame. Background from continuum processes and feed down 
from ip(2S) have not been subtracted. The shaded region has 
the expected shape for a contribution from B — > x^K. 

(1.53±g;|| ± 0.27) x 1CT 3 respectively. The non-zero \c2 
production is an indication that the factorization model 
does not give a complete picture for charmonium pro- 
duction in i3-meson decays. The momentum spectra in- 
clude a large low momentum component, indicating ei- 
ther multibody final states or final states with higher 
resonant K* production. 
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